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Metamaterials are effectively homogeneous materials that display extraordinary dispersion. Negative index 
metamaterials, zero index metamaterials and extremely anisotropic metamaterials are just a few examples. 
Instead of using locally resonating elements that may cause undesirable absorption, there are huge efforts to 
seek alternative routes to obtain these unusual properties. Here, we demonstrate an alternative approach for 
constructing metamaterials with extreme dispersion by simply coiling up space with curled channels. Such a 
geometric approach also has an advantage that the ratio between the wavelength and the lattice constant in 
achieving a negative or zero index can be changed in principle. It allows us to construct for the first time an 
acoustic metamaterial with conical dispersion, leading to a clear demonstration of negative refraction from 
an acoustic metamaterial with airborne sound. We also design and realize a double-negative metamaterial 
for microwaves under the same principle. 



Double negativity, zero index and extreme anisotropy are some of the unusual types of dispersion relations 
that can be realized using metamaterials, leading to applications including super-resolution imaging, 
invisibility cloaking and various transformation-based devices 122 . In particular, double-negative meta- 
materials, with separately tuned resonances in the two constitutive parameters, have been realized for electro- 
magnetic waves 2 5 , transmission-line circuits 6 , and acoustic waves 16 21 . These developments have stimulated the 
search for alternative routes to achieve negative refraction including chirality, band folding, time reversal, kinetic 
inductance and Mie scattering, which may provide a wider bandwidth, a lower loss, or a larger scalability to bulk 
metamaterials 23 31 . 

One of these alternative routes is to use a photonic or phononic crystal to achieve band folding; based on this 
approach, a band with negative phase velocity emerges that enables negative refraction 29,30 . Although photonic 
and phononic crystals lack a valid effective medium description and the wave manipulation is difficult to go 
beyond the diffraction limit, they have been proven to be a close approximation to achieve various extreme wave 
phenomena originally found in metamaterials. Apart from negative refraction, sharp focusing or imaging, zeroth- 
order Bragg photonic gap, invisibility cloaking and recently conical dispersion have been demonstrated in this 
way 31 " 37 . These wide-ranging demonstrations not only provide an alternative route to employing locally resonat- 
ing elements in conventional metamaterials but also an insight into the similarities between a photonic/phononic 
crystal and a metamaterial if we can reduce the band gap to low frequencies. Consequently, an effective medium 
can become meaningful and valid even at frequencies above the band gap. The extreme effective parameters at 
lower frequencies can allow us to manipulate waves in higher resolution and they can be further employed in the 
more general framework of transformation acoustics/optics in wave manipulations. 

In this case, if we wanted to avoid using locally resonating elements exhibiting high loss, a low-frequency band 
gap would normally require constitutive materials with high refractive indices, which may not be readily found in 
nature. For example, in the domain of acoustics for airborne sound, common solids have indices lower than one, 
while materials with a much lower sound speed than that of air (e.g. silicone rubber) usually accompany 
substantial absorption loss. It is also difficult to have high-index materials with low losses in the optical domain. 
As a result, the lattice constant binds to the working frequencies above the first band gap. Interestingly, in the case 
of metamaterials that rely on locally resonating units, there is also usually a binding between the electrical size of 
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the resonating units and the working frequency with negative indi- 
ces. We show that this binding can be relaxed so that the ratio of the 
wavelength to the lattice constant can be easily changed if we can 
yield a high refractive index by using a geometric route. 

In this article, we experimentally demonstrate a geometric 
approach for constructing metamaterials with extreme dispersion 
by coiling up space. By delaying the propagation phase using curled 
channels to mimic an array of high-index elements, an acoustic 
metamaterial with conical dispersion and associated double negativ- 
ity can be obtained at very low frequencies. As a result of our 
approach, we can demonstrate negative refraction from an acoustic 
metamaterial with airborne sound. Moreover, we also construct a 
double-negative metamaterial that works for microwaves under the 
same geometric principle. 

Results 

Design of space-coiling metamaterials. In the current scheme, the 
mentioned high refractive index is mimicked by coiling up the wave 
propagation space through curled channels with substantial phase 
delays (Fig. la). The refractive index n lD of such a one-dimensional 
element is the ratio of the actual elapsed phase to the corresponding 
elapsed phase in a straight channel of the same physical length. 
Suppose we join these elements into a two-dimensional (2D) array. 
The result is that band-folding occurs at a low frequency and that the 
array has an effective medium description with extreme indices 
around the band-folding frequency at the zone center 38 . The 
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Figure 1 | Principle in coiling up the acoustic and the electromagnetic 
space, (a) Waves confined along a curled channel have an equivalent 
refractive index n 1D ^s>l comparing to propagation along a straight line, 
(b) (photograph) A unit of the acoustic space-coiling metamaterial 
fabricated using aluminum with a = 2.33 cm. (c) Simulated pressure field 
in a single unit being excited by a plane wave propagating to the right at 
2.7 kHz, with red/blue color denoting high/low values, (d) (photograph) A 
unit of the electromagnetic space-coiling metamaterial by stacking copper- 
wire patterns printed on circuit boards, with a = 0.76 cm and foam spacers 
in the vertical direction, (e) Simulated H-field pattern at the copper wire 
layer being excited by a plane wave propagating to the right at 6.54 GHz. 
White lines represent the aluminum walls or the copper wires. 



wavelength ().) with these extreme indices can therefore be easily 
tuned; further, it can be much larger than the lattice constant (a) 
by simply geometrically increasing the number of turns of the curled 
channels. 

This space-coiling scheme applies to any situations in which waves 
can be guided without cutoff, including 2D and three-dimensional 
(3D) acoustic waves and 2D electromagnetic waves. Fig. lb shows a 
unit for a 2D acoustic space-coiling metamaterial (of a = 2.33 cm), 
fabricated by computer numerical control (CNC) milling on alumi- 
num. Fig. lc shows the full-wave simulation of an acoustic pressure 
field for a plane wave at 2.7 kHz (or a free-space wavelength of 
12.7 cm) impinging from the left. The channels (of width 1.5 mm 
and walls of width 0.8 mm, height is 1 cm) guide the acoustic waves 
to propagate in a curled fashion. From one corner to the center of the 
unit, the elapsed phase is found to be delayed by a factor of n 1D = 4.8, 
compared with the case when the two end points are connected by a 
straight line. In another case of electromagnetic waves of H-polar- 
ization propagating in 2D, Fig. Id shows a unit of the metamaterial 
(of a = 0.76 cm) designed for microwaves. The acoustic hard walls 
become the copper- wire (of width 0.15 mm and air-channel width of 
0.62 mm) layers printed on circuit boards (FR-4 PCB laminate). 
Identical patterns are then stacked layer-by-layer separated by foam 
spacers of 2.3 mm thickness. Full-wave simulation (performed at 
6.54 GHz) shows a similar wave confinement of the H-field in the 
curled channels (Fig. le). 

Conical dispersion and double negativity of space-coiling meta- 
materials. In the first step, we experimentally obtained the effective 
medium of space-coiling metamaterials. For an acoustic space- 
coiling metamaterial, we measured the complex transmission and 
reflection coefficients at normal incidence by inserting the meta- 
material between two waveguide tubes. The effective density (p), 
bulk modulus (B) and refractive index (n) were then extracted 
using the S-parameter retrieval method 3 '. The effective refractive 
index n from 2.7 to 3.7 kHz varied continuously from — 1 to 1 (cir- 
cles in Fig. 2a). The effective density p and reciprocal bulk modulus II 
B changed sign at around the same frequency (3.12 kHz; circles in 
Figs. 2b and 2c), corresponding to a normalized frequency Q 0 = mal 
2nc = l/«i£» in which the bands fold back to the zone center 38 . We 
also extracted the effective medium parameters from full-wave 
simulations with microstructures included using the same method 
(Fig. 2 solid lines), which showed very good agreement with the 
experimental results. Double negativity occurs below Q 0 > while the 
dispersion around Q 0 is conical in shape (upper inset of Fig. 2) 
with a refractive index that can be analytically expressed by 
n~\[2{nm — iKc/aja), giving rise to a nearly constant group index 
V2n 1D - 

While there have been previous demonstrations/proposals of con- 
ical (Dirac-cone) dispersion in the intermediate frequency regi- 
me 35 3y , our geometric route further relaxes the binding between X 
and a and makes possible a conical dispersive metamaterial. Indeed, 
the conical dispersion presented herein is highly isotropic and sym- 
metric in frequency. The group index of the dispersion found from 
the simulation (blue line, agreeing with the open circles obtained 
from the experiment in the upper inset of Fig. 2) deviates from the 
constant \/2nw by only 5% in any direction, even at a working 
frequency of a 15% increase/decrease from Q 0 . In fact, a tiny gap 
opens up at the cone center frequency because of the finite width of 
the air channel in constructing the unit cell 38 . However, the size of the 
gap is less than 1% of the cone center frequency. The analytic formula 
of the refractive index (a conical description) is thus an accurate 
representation of the dispersion, and agrees well with the experi- 
mental dispersion. 

Although our scheme is geometric in nature in order to avoid 
resonant absorption, there is still some practical loss. We thus experi- 
mentally extracted the imaginary part of the index (Im(n), circles in 
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Figure 2 | Effective medium of the space-coiling metamaterials. (a-d) The effective parameters of acoustic space-coiling metamaterial: the real part of 
effective index, density, reciprocal bulk modulus, and the imaginary part of effective index, (e-h) The effective parameters of electromagnetic space- 
coiling metamaterial: the real part of effective index, permittivity, permeability, and the imaginary part of effective index. Upper/lower inset: dispersion 
for the acoustic/electromagnetic metamaterial. The blue solid lines/open circles represent the simulation/experimental results. 



Fig. 2d). This part agrees with the full- wave simulation results (blue 
line in Fig. 2d) and is further confirmed by the smoother spectrum 
obtained from directly measuring the attenuation by cascading a 
variable number of unit cells (Fig. 2d, black solid line; see 
Supplementary Information). This absorption loss mainly comes 
from the viscosity of air within a thin air layer of skin depth 
(5 oc 1 / y^ai on the inner surfaces of the channels 40 . Moreover, it pro- 
vides a frequency dependence of Im(«)cc<5/flcc l/y/w/a, which 
agrees with both the experimental and the simulation results (blue 
line). The current acoustic sample has an experimental figure of 
merit (FOM = |Re(n)|/Im(n)) of 5.2 at 2.7 kHz with n = -1. If 
the size of the metamaterial is scaled up by s, the FOM can be further 
increased by ^/s times for a fixed value of n (by working at the same 
normalized frequency in our scheme). 

Next, we applied the space-coiling principle to electromagnetic 
waves. We measured the S-parameters by putting our electromag- 
netic (EM) space-coiling metamaterial (Fig. Id) between two stand- 
ard horn antennas for the H-polarization at normal incidence. Fig. 2e 
shows the experimentally extracted effective index (circles), together 
with the index obtained from the full- wave simulations (blue lines) 
with good agreement, varying from — 1 to 1 from 6.54 to 7.95 GHz, 
again indicating the occurrence of band folding at a very low fre- 
quency. Simultaneously negative e and (i were also found in the 
negative index regime (Figs. 2f and 2g). From the acoustic to the 
EM sample, the acoustic hard walls were converted into stratified 
layers of copper wires of the same pattern, separated by foam spacers 
of a subwavelength thickness 2.3 mm. The waves were then guided as 



TEM modes between two adjacent copper wires on the same layer. 
The H-polarization waves outside the metamaterial can then be 
coupled to the curled channels for propagation on each layer. 

Although the same principle induces double negativity/positivity 
regimes here again, there is a departure from the ideal conical dis- 
persion in the electromagnetic dispersion (lower inset of Fig. 2). A 
spectral gap now opens up from 7.1 to 7.4 GHz (a distinct phase of 
Re(n) ~ 0 in Fig. 2e or between the two zero-crossing frequencies for 
Re(s) and Re(|.i) in Figs. 2f and 2g). The foam spacer layers with finite 
thickness splits apart the effective electric and magnetic plasma fre- 
quencies. This gives us an additional way for tuning the conical 
dispersion with a variable gap size. By contrast, decreasing the thick- 
ness of the foam spacers would diminish the spectral gap (but this 
would also require more copper layers to construct the metamaterial 
as a compromise). In fact, the finite thickness of the foam spacer also 
introduces spatial dispersion into the flat band of electric plasma (e ~ 
0). Together with the mirror symmetry breaking of the wire pattern, 
the flat band can be excited in principle and additional sharp reso- 
nances are induced to the theoretical effective medium (Figs. 2f to 
2h). These sharp resonances are not significant here since they are 
smeared out by the tiny disorder of the sample, e.g. misalignment 
between the layers, in the experimental results. Thus, we have 
achieved a double-negative microwave metamaterial (isotropic in 
2D) with a small lattice constant around a < a/6, while the experi- 
mentally achieved FOM is around 10 (Figs. 2e and 2h at n ~ —1), 
including the additional scattering noise owing to the open-space 
environment (a fluctuating spectrum in the experimental Im(n) 
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results). One advantage of our geometric approach in coiling up 
space is that the ratio Ala can be easily changed in principle. XI a 
can be further increased when we employ a larger number of turns 
in the wire pattern. For example, if we fix the working frequency for 
a target refractive index, e.g. — 1, a larger number of turns results in 
a smaller lattice constant a and vice versa (see the Discussion sec- 
tion). This provides a potential method to make metamaterials work 
at the deep subwavelength limits, as natural materials do, which is 
important if we want to manipulate waves with deep subwavelength 
resolution. 

Demonstration of negative refraction with a slab. As the conical 
dispersion here is isotropic in 2D, the effective medium parameters 
experimentally established for the space-coiling metamaterials at 
normal incidence are generally valid for oblique incidence. By 
constructing the unit cells into a slab, the negative/positive indices 
below/ above the cone center frequency can be revealed by a negative/ 
positive beam shift across a slab. The upper inset of Fig. 3 shows the 
schematic of this beam shift experiment. A beam with an incident 
angle 0 t impinges to the slab and the green dashed line separates the 
two regions of negative and positive refractions. For acoustics, we 
construct and put a slab with 2 X 20 units (6.6 X 66 cm) in the 
middle of a 2D waveguide. An acoustic beam of 20 cm in width 
impinges the slab at 0, = 30°. The pressure field profile (20 cm 
away from the slab) is then scanned by a microphone (see Figs. 3a/ 
3b for the measured/simulated results) from 2.65 to 3.75 kHz. These 
two results actually agree very well in terms of beam positions and 



intensities. From the results, we can divide the frequencies into three 
regimes by two frequencies (the two yellow dotted lines) having n = 
±0.5, indicating the onset of total internal reflection. From 2.65 
to 2.89 kHz, the beam center appears on the left-hand side of 
the vertical green line. This indicates a negative index for our 
metamaterial. On the contrary, the beam center appears on the 
right-hand side for positive indices from 3.4 to 3.75 kHz. In 
particular, the beam intensity profiles at 2.7 and 3.7 kHz where n 
~ — 1 and 1 are shown in Fig. 3c, confirming the agreement between 
theory and experiment at both negative and positive refractions, (see 
Fig. S4 in Supplementary Information for the 2D pressure field 
profiles measured at these two frequencies.) In fact, the absorption 
inside the slab dampens the effect of multiple reflections, making the 
beam position generally follow Snell's law employing effective 
complex refractive indices (blue dashed line) while the intensity 
fluctuation at different frequencies persists (in both experiment 
and simulation) because of the impedance mismatch between air 
and metamaterial. The absorption also contributes to part of the 
amplitude drop (see Fig. S3 in Supplementary Information). From 
2.89 to 3.4 kHz, frustrated total internal reflection occurs. The beam 
tunnels through the barrier with a very small intensity except near 
the cone center. In this regime, the transmission coefficient can be 
approximated as t~(l — rf) exp( — |fc z |cf), where r s is the reflection 
coefficient of a single metamaterial surface, k z is the propagation 
constant in the direction normal to the slab and d is the thickness 
of the slab. Although the waves generally decay with an exponential 
factor, there is a resonance in the r s spectrum across the cone center 
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Figure 3 | Beam shift of a slab of space-coiling metamaterial. inset, schematic of beam shift experiment. A beam with incident angle 0 ; = 30° impinges to 
the slab. The transmitted beam position induced by negative indices locates to the left side of green dashed lines, and vice versa. Intensity profiles of output 
beam varying with frequency are measured, (a) Experimental and (b) simulation results for the acoustic slab of 2 X 20 units, (c) Selected normalized 
intensity profiles where n ~ — 1 at 2.7 kHz (blue) and n ~ 1 at 3.7 kHz (red), solid lines represent the simulation results, symbols represent the 
experimental results, (d) and (e) The electromagnetic cases corresponding to (a) and (b), the EM slab has 3 X 46 units X 68 copper layers, (f) Selected 
normalized intensity profiles of EM case where n ~ — 1 at 6.54 GHz (blue) and n ~ 1 at 7.95 GHz (red). Blue dashed lines in (a), (b), (d), (e) represent the 
beam center obtained by Snell's law with complex refractive indices. 
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frequency (owing to zero density there) and a small transmission 
around the cone center frequency results. 

Similarly, the beam-shift experiment for the electromagnetic 
space-coiling metamaterial, which spans 3 X 46 units in area with 
a height of 68 layers of copper wire patterns, is performed. A horn 
shines a beam (of width 9 cm) of H-polarization at 0, = 30°, while 
another horn scans the output beam profile at 16 cm away from the 
slab. The experimental result (Fig. 3d) agrees quite well with the 
simulation result with an effective medium extracted by the experi- 
mental process (Fig. 3e). There are some minor differences in the 
intensities (at ~ 7.7 GHz) because the experiment (for 2D wave 
propagation) was conducted in a 3D open space instead of a 2D 
waveguide as in the acoustic case. We can again divide the frequen- 
cies into three regimes, namely negative refraction (6.5 to 6.82 GHz), 
frustrated total internal reflection (6.82 to 7.62 GHz) and positive 
refraction (7.62 to 8 GHz). The experimental result still follows 
SnelFs law quite well. Specifically, the beam profile (blue) at 
6.54 GHz where n ~ — 1 highlights the negative refraction in 
Fig. 3f (theory and experiment) in contrast to the positive refraction 
at 7.95 GHz (n ~ + 1). In the tunneling regime, the resonance owing 
to the conical dispersion disappears because the waves exponentially 
decay more severely with a larger slab thickness in terms of wave- 
length in this case. 

Negative refraction from a prism. The isotropic effective medium is 
also valid for other geometries. We assembled acoustic unit cells into 
a 45° prism with a size of 4 X 4 units (Fig. 4a). A beam with width 
13.2 cm was then impinged normally on one surface and the 
refraction angle measured at a radius of 15 cm. From 2.7 to 
3.7 kHz, the direction of the experimentally measured refracted 
beam rotates in an anticlockwise direction, showing a continuous 
variation from a negative to a positive index (Fig. 4b). This is in good 
agreement with the simulation results shown in Fig. 4c, confirming 
the negative index of the sample as well as the validity of our 
experimentally extracted effective medium. 

Discussion 

It is interesting to note other approaches to obtaining a conical 
dispersion 35-37 in the intermediate frequency regime. Since our geo- 
metric approach reduces the first band gap to a very low frequency, it 
makes a conical dispersive metamaterial possible to construct. This 
also results in a conical dispersion that is highly isotropic and sym- 
metric in frequency. The group index fluctuates less than 5% even 
with a 15% increase/decrease in frequency from the cone center 
frequency. Actually, the ratio of the wavelength to lattice constant 
can be easily changed in principle in our geometric approach. In 
conventional metamaterials that rely on the local resonances of sub- 
wavelength structural units, the size of these structural units has to be 
a certain fraction of the wavelength to support resonances. By con- 
trast, II a here can be made either smaller or larger by changing the 



number of turns in coiling up the space. By using the microwave 
space-coiling metamaterial as an example, and by adding four more 
turns for the subwavelength channels in each unit (simulation in Fig. 
S8 in Supplementary Information), the ratio XI a (at a fixed frequency 
with n = —1) increases from 6 to around 8.1 (35% increase). 

In summary, we have demonstrated a geometric approach 
through coiling up space in order to achieve an acoustic metamaterial 
with the conical dispersion of constant group velocity and associated 
double negativity. The presented approach enables us the direct 
manifestation of negative refraction for isotropic acoustic metama- 
terials with air-borne sound. The same geometric principle is also 
valid for designing a double-negative microwave metamaterial where 
cutoff-free waveguiding can occur, while the ratio between wave- 
length and lattice constant can be easily changed in principle. 
These investigations will prove to be useful for designing future 
metamaterials with extreme dispersion. For example, a metamaterial 
with conical dispersion can be used for controlling slow wave pro- 
pagation in a broadband of frequencies with constant group velocity. 
By having a metamaterial with wavelength much larger than the 
lattice constant, we can manipulate waves in a finer resolution and 
the extreme effective parameters will be useful for employment in the 
framework of transformation acoustics or optics. 

Methods 

Acoustic wave measurement. We measure the S-parameters of the acoustic space- 
coiling metamaterial unit within a waveguide, with a rectangular cross-section of 
3.3 cm (width) X 1 cm (height), in order to extract the effective medium parameters 
in the time domain. A speaker is driven by a pulser-receiver (Olympus PR5800) on 
one side, while the transmitted or reflected signal is taken by an electret microphone, 
which is connected to an oscilloscope (Tektronix TDS1012B). Absorbers are set in the 
two open ends of the waveguide system to minimize reflections from the tube to the 
open space. The complex transmission coefficient is obtained by measuring the 
transmission signal through the metamaterial normalized to the case of air (without 
sample). In the reflection measurement, the source signal is subtracted by carrying out 
a control experiment without the sample while the complex reflection coefficient is 
obtained by measuring the reflection signal with reference to an aluminum block, 
which is considered to be a hard solid for sound waves propagating in air. For a 
pressure field-profile measurement, the device (see Supplementary Fig. SI for the case 
of a slab) is put into a 2D waveguide system using two parallel plates of tempered glass. 
A speaker fires a series of pulses to the device and a microphone is used to measure the 
pressure profile at different locations. A 2-axis automated positioning system from 
Velmex, Inc. is used to move the microphone in the waveguide through a pair of 
magnets with one side of waveguide sandwiched in between, (one magnet attaching to 
microphone within the waveguide and another attached to the positioning system) A 
similar setup is used to measure the negative/positive refraction from a prism 
constructed by the acoustic space-coiling metamaterial. 

Microwave measurement. The S-parameters are measured under the condition of 
normal incidence in an open space (Supplementary Fig. S2). The sample is located 
between a pair of aligned standard horn antennas (model Narda 642) and the 
measurements are taken for frequencies from 5.4 to 8.2 GHz. Linear H-polarization 
radiation is employed and a vector network analyzer (model E8361A PNA) is used to 
measure the S-parameters. S 2 i is normalized to that without the sample, while 
reflection Si i is normalized to the reflected signal from an aluminum plate with the 
same cross-section of the sample (18.4 X 49.7 cm). For beam shift measurement, we 
rotate the slab to make the incident angle 30' J . We then fix a horn as the source and 
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Figure 4 | Prism with acoustic space-coiling metamaterials. (a) A metamaterial prism with edge of 45 degrees inclination is impinged by a beam from 
the left-hand side. The (normalized) pressure profile of the output beam from the prism is measured at an arc of radius 15 cm from —60° to 60° 
(0" is defined along the green dashed lines), (b) and (c) The experimental and simulated pressure intensity profile, normalized to the peak intensity at each 
frequency. 
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horizontally move the other horn to scan the local power of the beam along the 
direction perpendicular to the previous propagation direction. We use a laser line 
level for alignment to keep the two horns aligned on a horizontal surface to fire 
H- polarization and to keep the sample interface vertical to the horizontal surface in 
order to allow us to investigate two-dimensional wave propagation. 

Simulation methods. We use the commercial software COMSOL Multiphysics for 
the full- wave simulations of both the acoustic and the EM space-coiling 
metamaterials. The simulations are performed including the microstructures with 
actual geometric sizes. The metals are modeled as ideal hard solids in acoustics and 
metals with finite conductivity in microwaves. The dielectric substrate for the PCBs 
(FR-4) is taken as £ - 3.8 + 0.03i. 
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